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(g) Verfahren und Vorrichtung zur oberfiachennahen Detektion von Stromdichteverteiiungen in einenn Uhtergrund 

(g) Die Erfindung betrifft ein Verfahren sowie eine Vorrich- 
tung zur oberfiachennahen Detektion von Stromdichte- 
verteiiungen in einem Untergrund, bei dem das in dem 
Untergrund verfahrensgemaB erzeugte Magnetfeld ge- 
messen und daraus die Stromdichteverteilung im Unter- 
grund bestimmt wird. Der Erfindung iiegt das technische 
Problem zugrunde, ein neues, nicht invasives Verfahren 
und eine entsprechende Vorrichtung anzugeben, die eine 
Bestimmung der Stromdichteverteilung mit groSerer Ge- 
nauigkeit ermoglichen. Erfindungsgemal^ wird bei dem 
Verfahren im Untergrund direkt ein Wechselstrom er- 
zeugt, dessen Magnetfeld vermessen wird. Dazu ist eine 
erfrndungsgema(le Vorrichtung vorgesehen. 
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Beschreibung 

Die Erfindung betrifft ein Verfahren sowie eine Vorrich- 
tung zur oberflachennahen Detekdon von Stromdichtever- 
teilungcn in einem Untergrund, bei dem das in dem Unter- 
grund verfahrensgemzlB erzeugte Magnetfeld gemessen 
wird und daraus die Stromdichteverteilung im Untergrund 
bestimmt wird. 

FCir das Verstandnis aller wesentlichen Fragen im Zusam- 
menhang mit Boden und Grundwasserleitem ist die Kennt- 
nis der rauniJichen Struktur des Untergrundes von entschei- 
dender Bedeutung. Diese Kenntnis wurde bisher vor allem 
durch invasive Verfahren, im wesentlichen durch Bohrun- 
gen, gewonnen. Diese Verfahren haben neben der Tatsache, 
dafi sie invasiv sind und die Verhaltnisse im Untergrund ver- 
andem konnen, den Nachteil, daB sie die raumliche Struktur 
nur punktuell an den Bohrungen wiedergeben. Die wichti- 
gen Zusammenhange und Verbindungen bleiben dabei un- 
entdeckt, die jedoch fur eine Gefahrdungsabschatzung, Vor- 
hcrsage und Beobachtung von Transponvorgangen im Un- 
tergrund notwendig sind. 

In j lingerer Zeit wurden zunehmend geophysikalische 
Methoden eingesetzt, die im wesentlichen nichi invasiv den 
Boden durchmustem. Hierbei liegen die Schwachen der 
Verfahren in der meist zu geringen raumlichen Auflosung, 
dem unklaren Zusammenhang der gemessenen physikali- 
schen GroBe mit dem eigentlichen interessierenden geologi- 
schen Parameter und der nicht eindeutigen Riickrechnung 
der mit den Sensoren aufgenommenen MeBwcrlc auf die zu- 
grundeliegenden Verhaltnisse im Untergrund. 

Ein solches, aus dem Stand der Technik bekanntes Ver- 
fahren wild als Geoelektrik bezeichnet, bei dem uber minde- 
stens zwei Elcktroden ein Gleichstrom in den Untergrund im 
Untersuchungsbereich eingeleiiet wird. Mil Elektroden und 
daran angeschiossenen Spannungsmefigeraten wird dann 
die resultierende elektrische Potential vertei lung an der Un- 
tergrundoberflache untersucht, woraus die elektrische Feld- 
verieilung im Untergrund bestimmt wird. 

Bei einem weiteren aus dem Stand der Technik bekannten 
Verfahren, das als Controlled Source Audio Magnetotelluric 
(CSAMT) bekannt ist, werden ein oder mehrcre elektrische 
Dipolc an den Boden angelegt und im Frequenzbereich von 
1 Hz bis 8 kHz betrieben. In einer groBen Entfemung (Gro- 
Benordnung etwa einige Kilometer) werden eine magne- 
tischc und eine elektrische Komponente der abgestrahlten 
elektromagnetischen Welle gemessen. Die vergleichsweise 
groBe Entfemung des Detektors ist erforderlich, um die ab- 
gestrahlte elektromagnetische Welle mathematisch mit der 
Femfeldnaherung beschreiben zu konnen. Daher werden 
mit dem genannten Verfahren geologische Strukturen in ei- 
ner groBen Tiefe und Ausdehnung gemessen. 

Weiterhin ist aus dem Stand der Technik ein induktives 
Verfahren der Geoelektromagnedk bekannt, bei dem die An- 
regung des Untergrundes mit Spulen an der Erdoberflache 
Oder mittels entfemter Sender bewirkt wird, indem im Un- 
tergrund als Reaktion auf das auBerste Magnetfeld im Unter- 
grund ein Magnetfeld induziert wird. Nachleilig bei diesem 
Verfahren ist, daB die Feldabschwachung mit steigendem 
Abstand in doppelter Weise auftritt, namlich einmal bei der 
Anregung des Untergrundes und dann bei der Messung 
durch den MeBsensor. 

Der vorUegenden Erfindung liegt daher das technische 
Problem zugrunde, ein neues, nicht invasives Verfahren und 
eine entsprechende Vorrichtung zur oberflachennahen De- 
tekdon von Stromdichteverteilungen in einem Untejrgrund 
anzugeben, wobei mit grofierer Genauigkeit die Stromdich- 
teverteilung im Untei;grund bestimmt werden kann. 

Die US 5,563,513 A offenbart ein Verfahren und eine 



Vorrichtung zur Widerstandsmessung in einem sedimentfor- 
migen Untergrund, wodurch anomale Widerstandsmuster 
nachgewiesen werden kormen. Die Vorrichtung wird auf ei- 
nem Untergrund angeordnet lind weist einen elektrischen 
5 Generator auf, der mit einem Transmitter ttiit geerdeten 
Elektroden verbunden ist. Zur Widerstandsmessung werden 
rechteckfbrmige Wellen von Strom mit groBer Amplitude 
und langen Perioden in den Untergrund eingeleitet, wodurch 
Magnetfeldveranderungen im Untergrund induziert werden, 
10 die an der Erdoberflache mit einem Magnetometer oder ei- 
ner Induktionsspuie gemessen werden konnen. Aus der Am- 
plitude und der Form des gemessenen Magnetfeldsignais, 
das als Funktion uber der Zeit aufgetragen wird, werden In- 
formationen iiber Widerstandsvariationen des Untergrund- 
15 aufbaus abgeleitet. 

Das zuvor aufgezeigte technische Problem wird erfin- 
dimgsgem^ durch ein Verfahren mit dem Merkmalen des 
Anspruches 1 gelost. 

GemaB dieses Verfahrens werden an dem Untersuchungs- 

20 bereich Elektroden mit dem Untergrund in Kontakt ge- 
bracht, durch die im Untersuchungsbereich ein elektrischer 
Wechselstrom hervorgerufen wird. Dazu wird iiber elektri- 
sche Zuleitungen eine elektrische Wechselspannung an den 
Elektroden angelegt. Durch unterschiedliche Leitfahigkci- 

25 ten in dem Untersuchungsbereich bildet sich eine entspre- 
chende Stromdichteverteilung im Untersuchungsbereich 
aus. Durch die Stromdichteverteilung wird ein Magnetfeld 
erzeugt, das mit Hilfe mindestens einen Magnetfeldscnsors 
an einer Mehrzahl von Orten gemessen wird. Die gemessene 

30 Magnetfeldstarkeverteilung wird dann mittels einer inver- 
sen Faltung mit dem Gesetz von Biot-Savat in eine Strom- 
dichteverteilung umgerechnet. 

In weiter bevorzugler Weise wird aus der Stromdichtever- 
teilung die Leilfahigkeitsverteilung im Untergrund berech- 

35 net. Dabei stellt die berechnete Leitfahigkeitsverteilung im 
Untergrund eine weitere ZielgroBe des Verfahrens dar, und 
die Interpretation der elektrischen Leitfahigkeitsverteilung 
im untersuchten Untexgrund kann dann abhangig von der je- 
weiligen Fragestellung in Bezug auf die hydraulische Leit- 

40 fShigkeit, die Wassersattigung oder die bevorzugten FlieB- 
wege abgeleitet werden. 

Die zuvor erw'ahnte inverse Faltung ist bei einem Pro- 
blem, das aufgrund der physikali schen Gegebenheiten als 
zwei dimensional anzusehen ist, eindeutig moglich. Fur die 

45 Inversion gibt es dabei eine Vielzahl von Moglichkeiten, 
beispielsweise auf der Basis von Finite-Elemente- oder Fi- 
nite-Differenzenmethoden, Fourieanalyse oder neuronalen 
Netzen. Dagegen ist eine eindeutige dreidimensionale In- 
version bisher nicht gelungen, jedoch werden aufgrund des 

50 erfindungsgemaBen Verfahrens Daten erzeugt, die eine 
merkliche Verbessenmg der Situation gegeniiber dem Stand 
der Technik herbeifuhren, so daB die Anzahl der fur eine 
dreidimensionale Inversion notwendigen sinnvollen Annah- 
men iiber die Stromverteilung im Untergrund verringert 

55 werden kann. 

Mit dem zuvor erlauteiten Verfahren lassen sich einerseits 
geologische und andererseits nicht geologische Untergriinde 
untersuchen. Als nicht geologische Untersuchungen konnen 
beispielsweise Werkstoffpriifungen zur Bestiirunung von 

60 Elissen oder Oxidadonsschaden in Metallen genannt wer- 
den. 

In einer ersten Ausgestaltung der vorliegenden Erfindimg 
wird bei der Berechnimg der Leitfahigkeitsverteilung aus 
der berechneten Stromdichteverteilung eine angenommene 
65 Verteilung des elektrischen Feides im Untergrund angenom- 
men. Dieses ist beispielsweise dann sinnvoU und moglich, 
wenn zwei flachige, plattenformige Elektroden im Unter- 
grund angeordnet sind, zwi schen denen sich eine nahezu ho- 
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mogene Potentialverteilung ausbildet. Daher kann in einem 
soichen Anwendungsfaii eine eleictrische Potentialvertei- 
lung wie bei einem Plattenkondensator in guter Naherung 
angenonunen werden. 

Bei weniger eindeudgen Potentialverteilungen kann das 5 
zuvor bescfaiiebene Verfahren dadurch unterstutet werden, 
dafi wahrend der Messung der MagnetfeldstSrkevcrteilung 
gleichzeitig die elektrische Potentialverteilung an der Ober- 
flache des Untergmndes gemessen wird. Somit konnen die 
gemessenen elektrischen Potentiale in die nahemngsweise lO 
Besdmmung der elektrischen Potentialverteilung mit be- 
riicksichtigt werden. 

Bei einer bevorzugten Ausgestaitung der vorliegenden 
Erfindung kann nait Hilfe einer Mehrelektrodenanordnung 
die raumliche Abhangigkeit der Leitfahigkeitsverteilung al- 15 
leine durcb mebrere Messungen der Stromdichteverteilun- 
gen fur zwei unabhangige Fliefirichtungen des Wechselstro- 
mcs, also ohne Messung oder Annahme bezuglich der Po- 
tentiale Oder der elektrischen Feider, bestimmt werden. We- 
gen der kleinen Frequenzen der Wechselspannung gilt nahe- 20 
rungsweise das Ohmsche Gesetz 

j = aE, 

wobei j der Stromdichievektor, a die (ortsabhangige) Leitfa- 25 
higkeit und E der Vektor des elektrischen Feldes ist. Das 
elektrische Fcld ist der Gradient des elektrischen Potentials 
0, also E = - V 0, folglich versch windet seine Rotation und 
es gilt: 

30 

0= V xE 
= V x(l/aj) 

= v(l/a)xj+ l/a(V xj) 

Damit ergibi sich die Vcktor-Gleichung: 35 

V xj + (a V 1/a) xj = 0 

Nach einer Messung der Stromdichteverteilung sind j und 

V x j bekannt, so da3 die letzte Gleichung als bestimmende 40 
Gleichung fiir 1/a angesehen werden kann. Das System ist 
noch unterbestimmt, weil wegen des Kreuzprodukts keine 
Aussage fiir die Komponenle von V (1/cr) parallel zum 
Stromdichtevektor gemacht werden kann. MiBt man aber 
(mindestens) zwei unabhangige Stromdichteverteilungen 45 
(z. B. durch Anderung der Elektrodenanordnung), so erhalt 
man ein iiberbesummtes Gleichungssystem, aus dem sich a 
absolut beslimmen lafit, wenn man die Leitfahigkeit an min- 
destens einer S telle, beispielsweise an den Zuleittmgen oder 

in einem wassergefiillten Hohlraum, kennt. 50 

Fiir die Losung des Gleichungssystems ist es giinsUg, 
wenn die beiden j -Feider annahemd senkrecht aufeinander 
stehen. Dies laBt sich z. B. dutch eine Vierelektrodenanord- 
nung emeichen, bei der die Verbindungslinien der in einer 
Messung verwendeten Elektrodenpaare senkrecht zueinan- 55 
der sind. 

Die Methode laBt sich insbesondere dann einsetzen, wenn 
Potentiaimessungen schwierig sind oder unzureichende In- 
formauonen liefem, z. B. wenn das zu unterscheidende Me- 
dium an der Oberflache unzureichend eiektrisch angekop- 60 
pelt ist Oder der zu untersuchende Bereich von einer sehr gut 
leitenden Schicht uberdeckt ist. 

Ein wesentlicher Vorteil des erfindungsgemafien Verfah- 
rens gegeniiber dem Verfahren der obengenannten Geoelek- 
trik besteht darin, dafi das Magnetfeld nicht durch leitMiige 65 
Bereiche im Untergrund abgeschirmt wird, so dafi auch Be- 
reiche des Unlergrundes unterhalb einer leitfahigen Schichi 
unlersucht werden konnen. Denn das Magnetfeld wird alien- 



falls durch ferromagnetische Einlagerungen beeiniiuBt. Die- 
ser Effekt kann jedoch - wie weiter unten erlautert wird - 
eliminien werden. Weiterhin liegt der Vorteil gegeniiber der 
Geoelektnk dahn, daB die Magnetfeldstarkeverteilung vek- 
toriell gemessen werden kann, wahrend die Messimg der 
Potentiale jeweils skaiarer Mefiwerte ergibt. Dadurch ist mit 
der Messung ein grofierer Infonnationsgehalt verbunden, 
wodurch auch die Inversion der gemessenen Magnetfeld- 
starkeverteiiung in eine Stromdichteverteilung detailreicher 
wird. 

Gegeniiber dem CSAMT- Verfahren, das oben beschrie- 
ben worden ist, liegt der Vorteil darin, daB in direkter Nahe 
zu dem Untersuchungsbereich das Magnetfeld vermessen 
wird, so daB der Untergrund erheblich detaillierter analysiert 
werden kann. 

Gegeniiber dem Verfahren der Geoelektromagnetik be- 
. steht der wesentliche Vorteil darin, daB der Strom mittels der 
Elektroden direkl in den Untergrund eingeleitet wird, ohne 
daB Spulen an der Untergrundoberflache oder entfemte Sen- 
der verwendel werden. Dadurch wird eine groBere Strom- 
dichte im Untergrund erteicht, so daB das dadurch erzeugte 
Magnetfeld genauer gemessen werden kann. 

In bevorzugter Weise wird die Magnetfeldstarke auBer- 
halb des Untergmndes gemessen, wobei insbesondere die 
Magnetfeldstarke in verschiedenen Hohen oberhalb des Un- 
tergnmdes gemessen wird. Dadurch besteht die Moglich- 
keit, die Veranderung des Magnetfeldes auch in vertikaler 
Richtung in verschiedenen Ebenen zu bestimmen, wodurch 
emeut die Dichte der aus der Messung abzuleitenden Infor- 
mationen gesteigert wird. Insbesondere gegeniiber der Geo- 
elektrik ergibt sich eine erhebliche Verbesserung, da die 
Geoeiektrik darauf beschrankt ist, die Feldverteilung an der 
Untergrundoberflache zu unlersuchen. Weiterhin ist es bei 
dem dargestellten Verfahren auch moglich, die Magnetfeld- 
starkeverteilung innerhalb des Untergmndes, vorzugsweise 
dabei in Bohrlochem, gemessen wird, Andererseils besteht 
die Mdglichkeit, bei Labomntersuchungen, bei denen der 
Untergrund sich innerhalb eines allseits zuganglichen Be- 
halters befindet, die Magnetfeldstarkeverteilimg seitli'ch ne- 
ben und/oder unteriialb der Laboranordnung zu messen. 

Die groBte Schwierigkeit bei der Messung von magne- 
tischen Feldem im Freien ist das standige Vorhandensein 
von Storfeldem. Die Herkunft der Slorf elder laBt sich in der 
folgenden Weise gliedem: 

- Das Erdmagnetfeld, das durch den EinfluB der Tro- 
posphare (Soimcnwind) zeitlichen Variationen unter- 
liegt. 

- Ferromagnetische Strukturen im Boden, die das Erd- 
magnetfeld vereerren. 

- Bewegte ferromagnetische Korper, die das Magnet- 
feld verzerren (2. B. vorbeifahrende Autos), 

- Antropogene Magnetfelder, etwa durch elektrische 
Maschinen, Uberlandleitungen, Radiowellen. 

- Induktion in bewegten eiektrisch leitenden Korpem. 

Die GroBe der Storfeider ist durchaus in der GroBenord- 
nung des zu messenden Signals aus dem Boden (einige nT), 
zum Teil, wie im Falle des Erdmagnetfeldes (50 pT), sogar 
einige Gr6Benordnungen grbBer. 

Die bevorzugte Auspragung des Verfahrens zur Kompen- 
sation dieser Storfeider liegt in der Verwendung von Lock- 
In-Technik. Der fiir die Untersuchung verwendete Strom 
wild mit einer festen Frequenz durch den Unto^uchimgsbe- 
reich geleitet. Die Messung des Magnetfeldes wird auf diese 
Frequenz eingeloggt, das heifit, es wird nur der Teil des Si- 
gnals vermessen, der dieselbe Frequenz wie der Untersu- 
chungsstrom aufweist. Auf diese Weise laBt sich das unter- 
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suchte Signal urn einige GrGBenordnungen vom Untergrund 
abheben. Ein typischer Wert fiir die Seitenbandunterdruk- 
kung (Dampfung des Signals bei halber/doppelter Frequenz) 
liegt bei 24 dB. 

Man kann davon ausgehen, dafi alle nicht durch den Un- 5 
tersuchiingsstrom hervorgerufenen Magnetfelder durch 
diese Methode herausgefiltert werden. 

Dabei kann die Lock-In-Technik zur Bestimmung von 
imaginaren Leitfahigkeitsfeldem herangezogen werden. 
Denn durch die Variation der Frequenz des verwendeten lO 
Stromas koimen die kapazitiv und induktiv hervoigenifenen 
imaginaren Leitfahigkeiten voneinander getxennt werden. 
Diese Trennung der imaginaren Widerstande ermoglichi 
Aussagen iiber die Porenstruktur und die innere Oberflache 
des Untergrundes. Der Effekt beruht dabei im wesentlichen 15 
auf Polarisation, einerseits fiir sehr gut leitende Inhomoge- 
nitaten durch die Polarisation des ionenleitenden Fluids und 
anderorseits fur schlecht leitende Bereiche durch die Polari- 
sation der Komer. Durch die Polarisation wird ein elektri- 
sches Feld aufgebaut, in dem Energie gespeichert ist, die 20 
dann zeitverzogert wieder abgegeben wird. Diese Zeitverzo- 
gerung laBt sich einerseits aus der Abklingkurve beim Ab- 
schalten bestimmen oder auch aus der Phasenverschiebung, 
wenn mit Wechselstrom gearbeitet wird. Letztendlich han- 
delt es sich hierbei um den kapazitiven Widcrstand als Teil 25 
des komplexen Widerstandes, fiir den ebenfalls die raumli- 
che Verleilung besummt werden kann. Die Interpretation 
dieser Information erstreckt sich vor allem auf das Auffin- 
den von Grenzfiachen unterschiedlicher Art der clcklrischen 
Leitfahigkeit (elektronisch, ionisch, durch Polarisation). Im 30 
Sinne geologischer oder transportrclevanter Parameter 
konnte man etwa auf Tone schlieBen, da besonders kleine 
Komer eine hohe Kapazitat hervomifen, analog zum Plat- 
tenkondensator, wo die Kapazitat reziprok vom Abstand der 
Platten abhangt. Das Auftreten elektronischcr Leitung in lo- 35 
kalisierten Bereichen, also das Vorhandensein von Metallen, 
oder auch das Auffinden schlecht leitender, polarisierbarer 
Fliissigkeiten ist ebenfalls moglich. 

Die gcmessene Magnetfeldstarkeverteilung kann aller- 
dings nicht ausschlieBlich auf die Stromvertcilung im Unter- 40 
suchungsbereich im Untergrund zuriickgefuhn werden. Vor 
allem zwei Magnetfelder unterschiedlichen Ursprungs miis- 
sen bei der Auswertung der MeBergebnisse beriicksichtigt 
werden. Zum einen wird durch den in den Zuleitungen in 
den Elektroden fliefienden Strom ein Magnetfeld erzeugt, 45 
und -zum anderen sind die sekundaren magnetischen Ant- 
worten, im wesentlichen durch Ferromagnetika, auf das er- 
zeugte Magnetfeld zu beachten. 

Die aus der Geophysik schon langer bekannte Methode, 
durch induktive Anregung von Wrbelstromen im Unter- 50 
grund die Verteilung der elektrischen Leitfahigkeit zu ermit- 
teln, kann in vorleilhafter Weise mit dem vorliegenden Ver- 
fahren kombiniert werden, um weiterhin die Aussagefahig- 
keit der gemessenen Daten zu erhohen. Denn die Aufnahme 
des Antwortsignals aus dem Untergrund kann mit dem glei- 55 
chen Sensoraufbau wie bei dem zuvor beschriebenen Ver- 
fahren durchgefiihrt werden. Dadurch lassen sich die beiden 
Verfahren leicht kombinieren, indem ledigiich beide Anre- 
gungsmoglichkeiten installiert werden. Die Reaktion des 
Untergrundes auf die verschiedenen Anregungen ist jewcils 60 
unterschiedlich, so daB durch die Kombination der Verfah- 
ren unabhangige Informationen gewonnen werden konnen. 

Durch das induktive Verfahren werden lokale Wirbel- 
strome im Untei^rund induziert. Diese bilden magnetische 
Dipole mit kurzer Dipollange (bezogen auf den Abstand der 65 
Sensoren). Fur ein Dipolfeld fallt die Feidstarke mit der drit- 
len Potenz in Abhangigkeit von der Entfemung ab. 

Bei der Messung mit direkt eingelegtem Strom ist die 



Entfemung des Sensors miiunter kleiner als die Lange eines 
Stromfadens (bzw. in der gleichen GroBenordnung). Man 
kann das magnetische Feld dann nahezu als zylindersymme- 
trisch ansehen. Dafiir fallt die Feidstarke mit der ersten Po- 
tenz in Abhangigkeit von der Entfemung ab. 

Wenn man also abwechselnd rait der beiden Anregungs- 
methoden miBt und sukzessive den vertikalen Abstand der 
Abtastung des Magnetfeldes andert, laBt sich aus dem unter- 
schiedlichen Intensitatsveriauf der Tiefenkontrast der Unter- 
suchung verbessem. 

Eine weilere Variation dieser Anordnung liegt darin, auch 
die anregende Induktionsspule in unterschiedlichen Hohen 
iiber dem Untersuchungsfeld anzubringen. Dabei ist einer- 
seits ein Modus denkbar, bei dem die Abtastung entspre- 
chend in der Hohe mitwandert oder aber unabhangig in ver- 
schiedenen oder einer konstanten Hohe durchgefiihrt wird. 

Bei einer weiteren bevorzugten Ausgestaltung der vorlie- 
genden Erfindung wird dem Untergrund in dem Untersu- 
chungsbereich eine Flussigkeit mit erhohter elektrischer 
Leitfahigkeit, insbesondere eine Salzldsung, zugefiihrt, um 
deren Transport und Verdiiimung zeitlich verfolgen zu kon- 
nen. Dieses geschieht vor alfein im Hinblick auf den StofF- 
transport im Untergmnd. In dieser Weise sollen bevorzugte 
Fliefiwege aufgespurt werden und ihre Bedeutung fiir den 
Transport in Boden und Grundwasserleitem untersucht wer- 
den. 

Generell kann das Verfahren zur Untersuchung der raum- 
lichen Struktur elektrischer Leitfahigkeiten von Objekten, 
die mil elektrischem Strom beaufschlagt werden konnen, 
eingesetzt werden. y 

- Ermittiung der Konzentration geloster Substanzen m ^ 
Boden und Grundwasserleitem aus der Leitfahigkeits- 
verteilung. 

- Beobachtung chemischer Prozesse durch Messung 
der zeitlichen und raumlichen Variationen der Leitfa- 
higkeits verteilung (z. B. Fallungsreaktionen). 

- Aufspuren und Untersuchen von elektrischen An- 
omalien (metallischer Korper, Erze, Leitungen, Blind- 
ganger, archaologischer Objekte, ...)• 

- Erkundung und Monitoring bevorzugter FlieB- imd 
Transportwege in porosen Medien (z. B. Boden, Aqui- 
fere, Kluftgesteine, Laboranordnungen). 

- Verlagerung von Agrochemikalien. j 

- Unterstiitzung des Katastrophenmanagements bei 
Induslrieunfallen . 

- Gefahrdungsabschatzung fiir Induslriestandorte. 

- Monitoring bei Labor- und Freilandversucfasfeldem 
zur Untersuchung von Transportvorgangen in porosen 
Medien. 

- Bestimmung von Partikel- oder koUoidalen Trans- 
port. 

- Aufspiiren und Untersuchen von besonders schlecht 
leitenden Fliissigkeiten in Boden und Grundwasserlei- 
tem (durch die Verdrangung des Wasserfilms auf den 
Komem), etwa: 

- Schadstofifahnen (Ole, PAKs) 

- Saniemng von Boden 

- Monitoring von Deponien auf Undichtigkeiten 

- Monitoring von Deichen auf Durch weichung. 

- Erkundung von wassertragenden oder wasserfuhren- 
den Schichten in ariden oder setniariden Gebieten. 

- Erkundung von Salz-SuBwasserfronten im Unter- 
grund, in FluBsystemen oder in Agrookosystemen. 

Das oben aufgezeigte technische Problem, wird erfin- 
dvmgsgemaB auch durch eine Vorrichtung nach Anspmch y€ f(\ 
zirr Durchfuhrung des Verfahrens zur oberflachennahen De- 
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tektion von Leitfahigkeitsverteilungen in einem Untergrund 
gelost. In den abhangigen Vorrichtungsanspruchen sind eine- 
Reihe von spezifischen Merkmalen aufgefiihit, die im ein- 
zelnen anhand von Ausfuhrungsbeispielen, die in der Zeich- 
nung dargestellt sind, eriautert werden. In der Zeichnung 
zeigt 

Fig- 1 ein erstes Ausfiihrungsbeispiel der wfindungsge- 
mafien Vorrichtung mit zwei piattenfonnigen Elektroden, 
die in einem voigegebenen Absiand zueinander angeordnet 
sind, in einer schematischen Darstellung, 

Fig. 2 ein zweites Ausfiihrungsbeispiel einer erfindungs- 
gemaBen Vorrichtung, bei der vier stabformige Elektroden 
im Untergrund angeordnet sind, in einer schematischen Dar- 
stellung, 

Fig. 3 ein drittes Ausfuhrungsbei spiel einer erfindungsge- 
mafien Vorrichtung, bei der die Elektroden in Bohriochem 
angeordnet sind, in einer schematischen Darstellung, und 
zeigen 

Fig. 4—7 weitere Ausfuhningsbeispiele mit unterschiedli- 
chen Elekirodenanordnungen und -ausgestaltungen. 

In Fig. 1 ist ein erstes Ausfiihrungsbeispiel einer erfin- 
dungsgemaBen Vorrichtung dargestellt, die einen Generator 
2 zum Erzeugen einer Wechselspannimg, eine Elektroden- 
anordnung mit zwei Elektroden 4, die plattenformig ausge- 
bildet sind und mit einem Untergrund 6 in Kontakt stehen, 
und Zuleitungen 8 aufweist, die die Elektroden 4 mit dem 
Generator 2 verbinden und somit einen Stromkreis herstel- 
len. Weiterhin weist die Vorrichtung eine Mehrzahl von Ma- 
gnctfeldsensorcn 10 zur Messung der Magnetfeldsiarkever- 
teilung inncrhalb eines Untersuchungsbereiches 12 auf, der 
in Fig. 1 mit einer gestrichelten Linie dargestellt ist. 

Bei dem in Fig. I dargestellten Ausfiihrungsbeispiel, das 
insbesondere fur Laboranordnungen geeignel ist, sind die 
Elektroden so angeordnet, daB sie wie bei einem Platten- 
Kondensator ein homogenes elektrisches Fcld erzeugen. Bei 
einer sol chen Anordnung laBt sich erstens sehr einfach von 
den Magnetfeldem auf die Strom verteilung schlieBen, da 
die Slromfaden in guter Naherung senkrecht zu den Elektro- 
den veriaufen und zweitens unmittelbar die Leitfahigkeits- 
verteilung abschatzen, da das elektrische Feld nahezu homo- 
gen ist. Mit numerischen Methoden konncn diese Abschat- 
zungen noch verbessert werden. 

Bei dem in Fig, 2 dargestellten zwei ten Ausfiihrungsbei- 
spiel, das insbesondere fur eine Anordnung im freien Ge- 
lande geeignet ist, sind zwei Paare von stabformigen Elek- 
troden 14 bzw. 16 im Untergrund 6 angeordnet, so daB sich 
eine gekreuzte Elektrodenanordnung ei]gibt. Der von dem 
Generator 2 erzeugte Wechselstrom fliefit bei zwei aufeinan- 
derfolgenden MeBreihen einmal zwischen den Elektroden 
14 und ein anderes Mai zwischen den Elektroden 16. Zu- 
satzlich lassen sich die Elektroden 14 imd 16 zu zwei unter- 
schiedlichen Quadropolen zusammenschalten, woraus wei- 
tere unabhangige Infonxiationen gewonnen werden konnen. 
Fiir jede Ausrichtung des zwischen den Elektroden 14 und/ 
Oder 16 fliefienden Wechselstromes wird die Magnetfeld- 
starkeverteiiung mit Hilfe der Magnetfeldsensoren 10 ge- 
messen. 

In Fig. 3 ist ein drittes Ausfiihrungsbeispiel einer erfin- 
dungsgemaBen Vorrichtung dargestellt, wobei - in perspek- 
tivischer Darstellung - insgesamt vier Bohrlocher IS darge- 
stellt sind, in denen jeweils Paare von Elektroden 20 ange- 
ordnet sind. Dadurch ist es mdglich, den Hauptteil des Hie- 
Benden Stromes durch den moglicherweise einige Meter un- 
terhaib der Uniergnindoberflache befindlichen Gnindwas- 
serleiter flieBen zu lassen. Daher kann die Stromeinleitung 
iiber Brunnen direkt in das Wasser erfolgen. Da das Wasser 
selbst eine um ein bis zwei GroBenordnungen groBere elek- 
trische Leitfahigkeit besitzt, als der wassergesattigte Boden 



typischerweise aufweist, muB der gesamte wassergefullte 
Bereich des Brunnens als KurzschluB fiir den Strom angese- 
hen werden. 

Wird die Stromeinleitung in den Untei^grund auf kurze 

5 Stiicke in jeweils groBer und geringer Hefe begrenzt, wie in 
Fig. 3 daiTgestellt ist, laBt sich auch die vertikale Korapo- 
nente der Stromfliefirichtung variieren, so daB die Datenba- 
sis der MeBwerte auch in dieser Richtung das Invertierungs- 
problem iiberbestimmt und eine vertikale DifFerenziening 

10 der Leitfahigkeits verteilung moglich wird. Erreicben laBt 
sich eine derartige Elektrodenanordnung, indem man die 
obere Elektrode an der nicht mit Grundwasser gesattigten 
Oberkante des Untergrundes anbringt und die untere Elek- 
trode innerhalb eines Brunnens so position! ert, daB die Elek- 

15 trode 20 nicht mit dem Wasser innerhalb des Brunnens in 
Kontakt steht, indem das Wasser oberhalb der Elektrode 
z. B. mit einem Packer herausgedriickt wird. Altemativ 
. kaim die Elektrode auch ohne Bohrung in den Untergrund 
gebracht (gerammt) werden, so daB sich oberhalb dieser 

20 kein wassergefuUter Hohlraum ausbildet. 

Eine weitere Ausgestaltung dieser Vorrichtung, die nicht 
in der Zeichnung dargestellt ist, basiert darauf, dafi sigmen- 
tiert aufgebaute Elektroden in den Untergrund eingebracht 
werden, wobei durch eine entsprechende Verschaltung der 

25 verschiedenen Elektrodenbereiche unterschiedliche Strom- 
flieBrichiungen realisiert werden konnen, 

Weiterhin besteht auch die Moglichkeit, durch unter- 
schiedliche Wahl der Geometrien der Elektrodenanordnun- 
gen den Fokus der durchgefiihrten Messung des Untersu- 

30 chungsbereiches zu verandem. Da sich die Stromdichtever- 
teilimg entlang der elektrischen Feldlinien im Untersu- 
chungsbereich ausbilden, ist in Bereichen hoher Feldlinien- 
dichte (Fokus) die Stromdichte besonders hoch und somit 
auch das hiervon ausgehende Magnetfeld. Dieses erhoht den 

35 Kontrast der Messung in diesen Bereichen hoher Strom- 
dichte. Durch eine variable Gestaltung der Geometrieen der 
Elektrodenanordnungen konnen demnach verschiedene Be- 
reiche des Untersuchungsbereiches detailliert untersucht 
werden. 

40 Die verwendeten Magnetfeldsensoren 10 konnen in un~ 
terschiedlicher Weise ausgebildet sein. Im wesentlichen sind 
dabei zu nennen: 

- SQUID (Superconducting Quantum Interference 
45 Device) 

- Induktionspulenmagnetomeier 

- MR (Magneto Resistance) 

- GMR (Giant Magneto Resistance) 

- Optisch gepulste Magnctometer/Gradiometer 
50 - Flux-Gate 

- Proton-Prazessionsmagnetometer 

Die MR/GMR Sensoren haben insbesondere den Vorteil, 
daB man sie zu groBen Arrays zusammenfiigen kann, wo- 

55 durch der MeBvorgang wesentlich vereinfacht wird. Die 
Squids zeichnen sich vor allem durch ihre extrem hohe 
Empfindlichkeit aus, sind durch die erforderUche Kiihlung 
allerdings schwer zu handhaben. 

Bei einer in den Fig. 1 und 2 nicht dargestellten Ausfuh- 

60 rung wird ein einzelner Magnetfeldsensor 10 verwendet, um 
das magnetische Feld oberhalb des Untersuchungsbereiches 
abzutasten. Wie in Fig. 2 dargestellt ist, besteht auch die 
Moglichkeit, eine Mehrzahl von Magnetfeldsensoren 10 
entlang einer Linie anzuoidnen, so daB ein eindimensionales 

65 Array von Magnetfeldsensoren 10 oreicht wird. Ein solches 
eindimensionales Array ISBt sich beispielsweise in Gestalt 
eines S tabes realisieren, entlang dessen die Magnetfeldsen- 
soren 10 angeordnet sind und wobei die Magnetfeldsenso- 
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ren 10 zusammen mit dem Stab uber den Untersuchungsbe- 
reich hin- und herbewegt werden. Dieses ist in Fig. 2 mit 
dem Doppelpfeil angedeutet. 

Weiterhin ist es auch mdglich, wie in Fig. 1 dargestellt ist, 
ein zweidimensionales Array aus Magnetfeldsensoren 10 zu 5 
bilden, das beispielsweise als Netz ausgefuhrt ist. Dieses 
zweidimensionale Array von Magnetfeldsensoren 10 wird 
auf dem Untersuchungsbereich 12 ausgebreitet, so daB an- 
schlieBend ohne Lageveranderung der Magnetfeldsensoren 
10 die Magnetfeldstarkeverteilung oberhalb des Untersu- lO 
chungsbereiches 12 mehrfach abgetastet werden kann. Das 
ist insbesondere dann von Voiteil, wenn man zeitabhangige 
Phanomene untersucht, indem eine Flussigkeit mit einer er- 
hohten elektrischen Leitfahigkeit, wie beispielsweise eine 
Saizlosung, auf den Untei^rund 6 aufgebracht wird und das 15 
zeitiiche Vertialten gemessen und analysiert wird. 

Fur ein genaues Vermessen der Magnetfeldstarkevertei- . 
lung ist die genaue Positionierung und Kontrolle der Aus- 
richtung des Magnetfeldsensors 10 bzw. der Magnetfeldsien- 
soren 10 erforderlich. Dabei kann in bevorzugter Weise der 20 
Sensor kardanisch aufgehangt sein und in Verbindung mit 
einer sich schnell drehenden Scheibe in einer Richtung fi- 
xiert werden. Zum anderen kann der Magnetfeldsensor 10 
mit einem Sender ausgestaliet sein, der von verschiedenen 
fest am Rand des Untersuchungsbereiches 12 installierten 25 
Empfangem angcpeilt wird, so dafi die Position des Magnet- 
feldsensors 10 immer bekannt ist. Die Anwendung weiterer, 
aus dem Stand der Tcchnik bekannler Verfahren zur Posiu- 
ons- und Lagebeslimmung konnen ebenso angewendet wer- 
den. 30 

Wie in den Fig. 1 und 2 dargestellt ist, sind die Zulcitun- 
gen 8 in besondercr Weise von dem Untersuchungsbereich 
12 wcggefuhn angcordnet. Dieses dient dazu, die von den in 
den Zuleitungen 8 flicBendcn Stromcn erzeugten Magnetfel- 
der zu kompcnsicrcn. Dazu werden die Zuleitungen 8 zu den 35 
Elektrodcn 4 bzw. 14 und 16 in Richtung des zu messenden 
Stromes vom Untersuchungsbereich 12 weggefuhrt und 
dann in rechien Winkeln in groBem Abstand zu einem 
Stromkreis verbunden. Fur cine derart kontrollierle Anord- 
nung der Zuleitungen 8 ist eine rechnerische Kompensation 40 
der Magnetfelder nach der Messung moglich, wenn wah- 
rend der Messung der in den Zuleitungen 8 flieBende Strom 
gemessen wird. Hierzu muB zudem die induktiv angeregte 
Magnetfeldantwort beriicksichtigt werden. Dazu sind in be- 
vorzugter Weise die Zuleitungen 8 zumindest teilweise als 45 
Koaxialkabel 9 ausgebildet, so daB die erzeugten Magnetfel- 
der zumindest teilweise abgeschirmt werden. Dadurch laSt 
sich der Standort des Generators beliebig, auch innerhalb 
des Untersuchungsbereiches wahlen, s. Fig. 4. Weiterhin ist 
auch eine Abschirmung aus p-Melall moglich. 50 

Ebenso ist eine magnetische Abschirmung der Zuleitun- 
gen 8 mittels eines Supraleiters moglich, jcdoch wird da- 
durch der technische Aufwand insbesondere bei Freilandun- 
tersuchungen vergroBert. 

Weitere Ausgestaltungen der Zuleitungen sind in den Fig. 55 
5 bis 7 dargestellt. 

Fig. 5 zeigt ein Ausfiihrungsbeispiel, bei dem die Zulei- 
tung zu den Eleku-oden uber Koaxialkabel 9 und ein Gitter- 
netz 22 vorgenommen wird. Durch das Koaxialkabel 9 flie- 
Ben gleiche Strome in entgegengesetzter Richtung, so daB 60 
sich um das Kabel 9 die resultierenden Magnetfelder kom- 
pensieren. Beim Einsatz eines Gittemetzes 22 sind die 
Stromdichten durch die Zuleitung und dadurch die resultie- 
renden Storfelder kleiner als bei einem einzelnen dunnen 
Leiter geringer. AuBerdem laBt sich das Gittemetz 22 so 65 
auslegen, daB durch das Zuleitungsmagnetfeld die Feldver- 
teilung eines ideal homogenen Unteigrundes gerade kom- 
pensiert wiirde. Im realistischen Fall eines heterogenen Un- 
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tergrundes werden also die Kontraste hervorgehoben. Statt 
eines Gittemetzes kann auch eine kaschiene Leiterplatte 
eingesetzt werden. 

Fig. 6 zeigt ein Ausfiihrungsbeispiel, bei dem die beiden 
Elektrodenstrome zusammen uber parallele Kabel oder ein 
Koaxialkabel 9 transportiert werden, so daB sich kein Storsi- 
gnal ergibt. Daher ist es moglich, eine Eiektrode 16 mit dem 
Generator 2 zu verbinden, ohne Storsignale zu erzeugen. 
Das Kabel muB aber anschlieBend iiber die Leitung 8 mit der 
zweiten Eiektrode 16 verbunden werden, wobei hier nur 
eine Richtung ein Strom flieBt, und damit ein storendes Ma- 
gnetfeld erzeugt. 

Bei dem in Fig- 7 dargestellten Ausfiihrungsbeispiel wird 
auch dieses Magnetfeld korapensiert, indem parallel zur Zu- 
leitung 8 eine weitere Leitung 24 verlegt wird, uber die ein 
zweiter Generator 26 mit einem gegenphasigen Wechsel- 
strom einen Kondensator 28 am Ort der zweiten Elekttode 
18 auf- bzw. entiadt. Die resultierenden Felder der Zuleitun- 
gen 8 und 24 heben sich dadurch gegenseitig auf. Damit der 
Kompensationsstrom durch die Zuleitung 24 standig mit 
entgegengesetzter Phase, aber gleicher Amplitude flieBt wie 
durch Leitung 8, wird der Generator 26 uber ein StrommeB- 
gerat 30 mit dem Signal aus der Zuleitung 8 getriggert. Um 
eine moglichst vollstandige Kompensation zu erreichen, 
werden die Zuleitungen 8 und 24 in einem Koaxialkabel zu- 
sammengefaBt. 

Sind weiterhin die Magnetfeldsensoren 10 in einem zwei- 
dimensionalen Array netzartig angeordnet, so kann das Netz 
gleichzeitig zur Messung des Stromes, der das Ncu an den 
jeweiligen Stellcn durchflieBt, verwendct werden. Dadurch 
wird eine sehr genaue Kompensation der Zuleitungsfelder 
moglich. An das zuvor beschriebene Netz laBt sich zudem 
eine Vorrichtung zur Messung des elektrischen Potentials, 
das durch den StromfiuB im Untergrund 6 an der Gbcrflache 
des Untersuchungsbereiches 12 hervorgerufen wird, anbrin- 
gen. 

Bereits oben wurde darauf hingewiesen, daB das von dem 
zugeleiteten Wechselstrom im Untergrund hervorgcrufene 
Magnetfeld zudem durch moglicherweise vorhandenc ferro- 
magnetische Materialien beeinfluBt wird. Denn es konnen 
Suszeptibilitaten in der GroBenordnung von 10^ in reinen 
Ferromagneten auftreten. Daher wird das magnetische Si- 
gnal der Stromverteilung also von dem Magnetfeld lokali- 
sierter Dipole uberlagert, wenn ferromagnetische Materia- 
lien im Untergrund enthalten sind. Fur eine Kompensation, 
bzw. eine Untersuchung, der lokalisierten ferromagneti- 
schen Einlagerungen lassen sich im wesentlichen drci Kon- 
zepte angeben: 

1. Induktive Anregung 

Durch eine Leiterschleife, die auf dem Erdboden um 
das Untersuchungsgebiet gelegt wird, wird ein Ma- 
gnetfeld in den Untergrund eingekoppelt, das von den 
ferromagnetischen Bereichen verstarkt wird. Die Ma- 
gnetfeidverteilung wird in der gleichen Weise wie bei 
der eigentlichcn Messung gemessen. Die Lock-In- 
Technik kann auch hier angewendet werden. Dabei 
sind niedrige Frequenzen im Bereich kleiner als 1 Hz 
zu wahlen, um die Anregung von Wirbelstromen im 
Untergrund gering zu halten. Dadurch, daB praktisch 
keine Strome flieBen, wird auf diese Weise die Vertei- 
lung feiTomagnetischer Bereiche im Untergrund be- 
sdmmt, die anschlieBend bei der Messung mit Wech- 
selstrom herausgerechnet werden kaim. 
Bei der einfachsten Beschaltung (Zuleitungen zu den 
Stromeinleitungselektroden) auf dem Feld, ist durch 
die Zuleitungen bereits eine solche Leiterschleife gege- 
ben (s. Fig. 1 und 2). In diesem Fall gefat man so vor. 
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dafi mit zunachst oberirdisch kurzgeschlossenen 
Stromleiteiektroden die induktive Messung durchge- 
fiihn wird, um danach die Stronunessung durchzufiih- 
ren. 

2, Da es sich bei den ferromagnetischen Bereichen um 
induktive Leitfahigkeiten bzw. Widerstande handelt, 
bewirken diese eine frequenzabhangige Phasenver- 
schiebung des MeBsignals an dieser Stelle. Die Fre- 
quenzabhangigkeit ist umgekehrt zu derjenigen des ka- 
pazitiv bedingten imaginaren Wdersiandes, so dafi 
diese beiden moglichen Quellen von Phasenverschie- 
bungen durch Frequenzvariation getrennt werden kon- 
nen (kapazitiver Widerstand sinkt mil der Frequenz: 
1/iwc, induktiver Widerstand steigt mit der Frequenz 
iwL). 

3. Eine dritte Moglichkeit der Kompensation erwachst 
aus der Tatsache, dafi es sich bei den Einlagerungen im . 
wesentlichen um lokalisierte Bereiche handelt, die ent- 
sprechend ein Dipolfeld mit kurzer Dipollange erzfeu- 
gen, Verglichen mit der durch den Stromflufi hervoige- 20 
rufenen Feldverteilung ist also mit viel hdheren Feld- 
kriimmimgen zu rechnen. Durch die rechnerische An- 
passung von Dipolfeldem an entsprechende Bereiche 
der gemessenen Werte konnen die Einfliisse ferroma- 
gnetischer Bereiche dann aus dem Signal der Magnet- 25 
feldverteilung herausgerechnet werden, bevor die In- 
version durchgefiihrt wird. 



Patenianspriiche 

1. Verfahren zur obcrflachennahen Detektion von 
Stromdichteveneilungcn in eincm Untergrund, 

- bei dem mindestens zwei Elektroden von min- 
destens einer Elektrodenanordnung mit dem Un- 
tcrgrund in Kontaki gebracht werden, 

- bei dem iiber elektrische Zuleitungen eine elek- 
trische Wechselspannung an den Elektroden angc- 
legt wird, 

- bei dem sich im Untergrund aufgrund der Leit- 
ffr^t.iAC^Lc^^S^l:;'<l^^ ^^"^ Stromdichteverteilung 

enct.*.^^ oausbreitet, 

_ - bei dem die Magnetfeldstarkeverteiiung, die 
ri25i^4^^<r^^ -durch die im Untergrund flieBende Stromdichte- 
'fc^le* Oy^jr verteilung erzeugt wird, mit Hilfe mindestens ei- 
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^Uii /pu ^ Xt- Magnetfeldsensors an einer Mehrzahl von Or- 
aJ'^^^i r**^ *^>i ten innerhalb eines Untersuchungsbereiches ge- 
pd^alU V&v<:'<1-n,essen wini, und 

K\f <^ ~ bei dem die gemessene Magnetfeldslarkevertei- 

o J \\xvi% mittels einer inversen Faltung mit dem Ge- 

setz von Biot-Savart in eine Stromdichtevertei- 
lung umgerechnel wird. 
_ 2. Verfahren nach Anspruch 1 , bei dem aus der bcrech- 
5T) V -^rW io^^i^ieten Stromdichteverteilung die Leitfahigkeitsvertei- 
V X ^ 4ung im Untergrund berechnet wird. 

aZ^.|o*imt^'^ 3_ Verfahren nach Anspruch 2, bei dem die Leitfahig- 
keitsverteilimg unter Zuhilfenahme einer angenomme- 
nen Verteilung des elektrischen Feldes im Untergrund 
berechnet wird. 

4. Verfahren nach Anspruch 3, bei dem wahrend der 
Messung der Magnctfeidstarkeverteilung gleichzeitig 
die elektrische Potentialverteilung an der Oberflache 
des Untergrundes gemessen wird. 

5. Verfahren nach Anspruch 2, bei dem fiir mindestens 
zwei verschiedene Elektrodenanordnvmgen die Mes- 
sung der Verteilung der Magnetfeldstarke durchgefiihrt 
wird, wobei die Hauptausrichtungen der Stromdichte- 
verteilungen fiir jeweils zwei Elektrodenanordnungen 
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unter einem vorgegebenen Winkel, vorzugsweise von 
90°, zueinandor verlaufen, und daB aus der berechneten 
Stromdichteverteilung dirckt die Leitfahigkeitsvertei- - 
lung berechnet wird. 

6. Verfahren nach einem der Anspriiche 1 bis 5, bei 
dem die Magnetfeldstarke aufierhalb des Untergrundes 
gemessen wird. 

7. Verfahren nach Anspruch 6, bei dem die Magnet- 
feldstarke in verschiedenen Hohen oberfaalb des Unter- 
grundes gemessen wird. 

8. Verfahren nach einem der Anspriiche 1 bis 7, bei 
dem die Magnetfeldstarkeverteilimg innerhalb des Un- 
tergrundes, vorzugsweise in Bohrlochem, gemessen 
wird. 

9. Verfahren nach einem der Anspriiche 1 bis 8, bei 
dem als Untergrund eine Laboranordnung verwendet 
wird und die Magnctfeidstarkeverteilung iiber und/ 
Oder seitlich neben und/oder unterhalb der Laboranord- 
nung gemessen wird. 

10. Verfahren nach einem der Anspriiche 1 bis 8, bei 
dem als Untergrund ein zu uberpriifender Werkstoff 
verwendet wird: 

11. Verfahren nach einem der Anspriiche 1 bis 10, bei 
dem mittels Lock-In-Technik, die mit der Frequenz der 
Wechselspannung ausgefuhrt wird, die Magnetfeld- 
starke gemessen wird, und insbesonderc bei dem die 
raumliche Verteilung der Phasenlage zwischen ausge- 
pragter Spannung und gcmessenem Magnetfeld ausge- 
wcrtet und ein raumlichcs Bild der Polarisierbarkcit 
des Untergrundes ermittelt wird. 

12. Verfahren nach Anspruch 1 1 , bei dem die Messung 
der Magnetfeldstarke mit der Lock-In-Tcchnik zu un- 
tcrschicdlichen Frequenzcn der angel cgten Wechsel- 
spannung durchgcfuhrt wird. 

13. Verfahren nach einem der Anspriiche 1, bis 12, bei 
dem die durch den in den elektrischen Zuleitungen flic- 
Bendcn Strom hervorgcrufenen Storfelder bcstimmt 
und von den gemessenen Magnetfcldstaxken abgezo- 
gen werden. 

14. Verfahren nach einem der Anspriiche 1 bis 13, bei 
dem die durch ferromagnctischc Elemente im Unter- 
grund induzierten St5rfeldcr bcstimmt und von den ge- 
messenen Magnelfcldstarken abgezogen werden. 

15. Verfahren nach einem der Anspriiche 1 bis 14, bei 
dem mit Hilfe einer das Untersuchungsgebiet zumin- 
dest leilweise umfassenden Lei terse hleife ein Magnet- 
feld im Untergrund induziert wird und bei dem bei ver- 
schiedenen Abstanden der Leiterschleife und/oder des 
mindestens einen Magnetfeldsensors vom Untergrund 
das induzierte Antwortsignal gemessen wird. 

16. Verfahren nach einem der Anspriiche 1 bis 14, bei 
dem eine Flussigkeit mit erhohter elektrischer Leilfa- 
higkeit, insbesonderc eine Salzlosung, dem Untergrund 
zugefUhrt wird. 

17. Verfahren nach Anspruch 16, bei dem zu unter- 
schiedlichen Zeitpunkten vor und/oder nach dem Ein- 
leiten der Flussigkeit die Leitfahigkeitsverteilung im 
Untergrund bcstimmt wird und die Differenzen der ver- 
schiedenen Leitfahigkeitsverteilungen berechnet wer- 
den, 

18. Verfahren nach Anspruch 17, bei dem aus der Er- 
hGhung der Leitfahigkcit die Konzentration eines Elek- 
trolyten im Untergnmd beslimmt wird. 

19. Vorrichtung zur Durchfuhrung des Vcrfahrens zur 
obcrflachennahen Detektion von Stromdichteverteilun- 
gen in einem Untergrund, insbesonderc nach einem der 
Anspriiche 1 bis 18, 

- mit einem Generator (2) zum Erzeugen einer 
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Wechselspaimung, ^a) ^.i^i^z^Jti^ enf^f- /iSoc^^J^^Zi^i^*-^^ 

- mit mindestens einer Elektrodenanordnung mit ' ^ ' ^ / 
mindestens zwei Elektroden (4, 14, 16, 20), die ^) U-oc^c^. ^^^tMV^ 
mit dem Untergnind (6) in Kontakt stehen, 

- mit Zuleitungen (8), die die Elektroden (4, 14, 5 lit /? '^ i \ " ^ \ 
16, 20) mit dem Generator (2) verbinden, und ^m^K^ dj ^^t, Oj 

- mit mindestens einem Magnetfeldsensor (10) * ^ r 

zur Messung der Magnetfeldstarkeverteilimg in- d.^ T' I- S^u^o^ 

nerhalb eines Untersuchungsbereiches (12) a — ii^ ^"H y^ViC 

dadurch gekennzeichnet, daB ein Lock-In^Verstarker 10 fo(^ \JJi:,'Ac>^i^ h^^^ a n\ 

zur Auswertung des MeBsignals vorgesehen ist. 

20. Vorrichtung nach Anspruch 19, dadurch gekenn- 
zeichnet, daB die Zuleitungen (8) zu den Elektroden (4, 
14, 16, 20) in der Richtung des zu messenden Stromes 
vom Untersuchungsbereich (12) weggefuhrt und mit 15 
einem vorgegebenen Abstand zu dem Untersuchungs- 
bereich (12) in vorgegebenen Winkeln, vorzugsweise . , 
in rechten Winkeln, zu einem Stromkreis mit dem Ge- 
nerator (2) verbunden sind. 

21. Vorrichtung nach Anspruch 19 oder 20, dadurch 20 
gekennzeichnet, daB die Zuleitungen (8) zumindest 
teilweise als Koaxialkabel (9) ausgefuhrt sind, 

22. Vorrichtung nach einem der AnsprUche 19 bis 21, 
dadurch gekennzeichnet, daB die Elektroden (14, 16) 
Liber ein Gittemetz (22) oder eine kaschierte Leiter- 25 
plattc mil den Zuleitungen (8, 9) verbunden sind. 

23. Vorrichtung nach einem der Anspruche 19 bis 22, 
dadurch gekennzeichnet, daB die Elektrodenanordnung 
zwei flachigc, plaitenformigc Elektroden (4) in einem 
vorgegebenen Abstand, vorzugsweise parallel zueinan- 30 
der ausgerichlet, aufweisl. 

24. Vorrichtung nach einem der Anspruche 19 bis 22, 
dadurch gekennzeichnet, daB die Elektrodenanordnung 
zwei stabformige Elektroden (14, 16) mit einem vorge- 
gebenen Abstand aufwcist. 35 

25. Vorrichtung nach einem der Anspruche 19 bis 22, , 
dadurch gekennzeichnet, daB die Elektroden (20) in 
Bohrlochem (18) in den Untergrund (6) eingebrachi 
sind. 

26. Vorrichtung nach Anspruch 25, dadurch gekenn- 40 
zeichnet, daB die Elektroden (20) in unterschiedlichen 
Hohen innerhalb der Bohrlocher (18) mit dem Unter- 
gnind (6) in Kontakt stehen. 

27. Vorrichtung nach einem der Anspruche 19 bis 26, 

dadurch gekennzeichnet, daB ein einzelner Magnet- 45 - 
feldsensor (10) zu einer Mehrzahl von MeBstellen be- 
wegbar angeordnet ist. 

28. Vorrichtung nach einem der Anspruche 19 bis 26, 
dadurch gekennzeichnet, daB eine Mehrzahl von Ma- 
gnetfeldsensoren (10) entlang einer Linie angeordnet 50 
sind und die Magnetfeldstarke an einer Mehrzahl von 
Orten entlang einer Linie messen. 

29. Vorrichtung nach einem der Anspruche 19 bis 26, 
dadurch gekennzeichnet, daB eine Mehrzahl von Ma- 
gnetfeldsensoren (10) entlang einer Flache verteilt an- 55 
geordnet sind und die Magnetfeldstarke an einer Mehr- 
zahl von Orten auf der Flache messen. 
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1- Introduction. There are a variety of medical problems for which it would be 
useful to know the time-varying distribution of electrical properties inside the body. 

By "electrical properties," we mean speci cally the electric conductivity and per- 
mittivity. The electric conductivity is a measure of the ease with which a material 
conducts electricity; the electric permittivity is a measure of how readily the charges 
within a material separate under an imposed electric eld. High- conductivity ma- 
terials allow the passage of both direct and alternating currents; high- permittivity 
materials allow the passage of only alternating currents. Both of these properties are 
of interest in medical applications, because di erent tissues have di erent conductiv- 
ities and permittivities. 

One medical problem for which knowledge of internal electrical properties would 
be useful is the detection of pulmonary emboli, or blood clots in the lungs. The devel- 
opment of pulmonary emboli is a common, and often serious, complication of surgery. 
Unfortunately, at present the diagnosis is rather involved, requiring inhalation of ra- 
dioactive gas in order to determine the ventilated lung region. This is followed by 
injection of a radio-opaque dye or a dissolved radioactive substance into a vein to 
make an image of the blood circulation. The image of the circulation in the lung is 
compared with the image of the ventilated region; areas that are ventilated but not 
perfused by blood indicate the presence of emboli. 

However, another way to determine the location of gas ajid blood within the body 
would be to map the internal electric conductivity and permittivity. These electrical 
properties are very di erent for air, tissue, and blood; moreover, they vary on di erent 
time scales. Thus a time- varying map of the electrical properties should show lung 
regions that are ventilated but not perfused by blood. 

Received by the editors January 30, 1998; accepted for publication (in revised form) May 21, 
1998; published electronically January 22, 1999. This work was pauiiially supported by the National 
Institutes of Health, the National Science Foundation, and the O ce of Naval Research. 
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Furthermore, determining the presence of pulmonary emboh from a map of the 
electrical properties would have a number of advantages over present techniques. It 
would require no exposure to X-rays or radioactive material. It could be done at the 
bedside, with a relatively small and inexpensive electrical system. 

Information about the internal electrical properties of a body could have many 
other medical uses. Such information could potentially be used for the following [H]: 
monitoring for lung problems such as accumulating uid or a collapsed lung, nonin- 
vasive monitoring of heart function and blood ow, monitoring for internal bleeding, 
screening for breast cancer, studying emptying of the stomach, studying pelvic uid 
accumulation as a possible cause of pelvic pain, quantifying severity of premenstrual 
syndrome by determing the amount of intracellular vs. ^extracellular uid, determining 
the boundary between dead and living tissue, measuring local internal temperature 
increases associated with hyperthermia treatments, and improving electrocardiograms 
[CF] and electroencephalograms. 

A variety of nonclinical applications of electrical impedance tomography is also 
possible. These include imaging multiphase uid ow [XHHB, WB], determining the 
location of mineral deposits in the earth [DL, P, SSS], tracing the spread of contam- 
inants in the earth [RDLOC, RDBLR, DR], nondestructive evaluation of machine 
parts [ESIC], and control of industrial processes such as curing and cooking. 

In order to map the electric conductivity and permittivity inside the body, our 
group at Rensselaer has designed and built an electronic system that applies currents 
through electrodes attached to the surface of the body and measures the resulting 
voltages. This system uses these electrical measurements to reconstruct and display 
approximate pictures of the electric conductivity and permittivity inside the body. 
This process is called electrical impedance tomography (EIT). The term "impedance" 
comes from circuit theory: it is the ratio of the voltage across a circuit element to the 
current through the element. 

Below we describe the mathematical model for EIT. We use this model to describe 
some of the theory that gave rise to the design of the Rensselaier system, which we call 
the adaptive current tomography (ACT) system. Then we survey some reconstruction 
algorithms, focusing on the one that was used to make the images accompanying this 
paper. 

2. The Mathematical Model. The electric potential u in the body is governed 
by the equation 

(2.1) V (x,u;)Vu=0. 

Here x is a point in , u is the electric potential or voltage, and the admittivity is 
given by (x,u) = {x,oj) + iu (x^oj), where is the electric conductivity, is the 
electric permittivity, and cj is the angulax frequency of the applied current. It is the 
admittivity of a block of homogeneous material that is proportional to the reciprocal of 
its impedance. Appendix 1 shows how (2.1) can be obtained from Maxwell's equations. 

In practice, we apply currents to electrodes on the surface d of the body. These 
currents produce a current density on the surface whose inward pointing normal com- 
ponent is denoted by jf. Thus 

(2.2) ^ = j on d . 

One possible model for EIT is (2.1) and (2.2), together with the conservation of 
charge condition j = 0 and the condition u = 0, which amounts to choosing 
a "ground" or reference voltage. This model is the commonly used continuum model. 
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For the continuum model, we de ne the operator R by Rj = v, where v denotes 
the restriction of u to the boundary. 

Unfortunately, the continuum model is a poor model for real experiments [CING], 
because we do not know the current density j. In practice, we know only the currents 
that are sent down wires attached to discrete electrodes, which in turn are attached 
to the body. One might approximate the unknown current density as a constant 
over each electrode (the gap model) ^ but this model also turns out to be inadequate 
[CING]. 

We need to account for two main e ects: the discreteness of the electrodes, and 
the extra conductive material (the electrodes themselves) we have added. We can do 
this as follows [CING, SCI]. 

The integral of the current density over the electrode is equal to the total current 
that ows to that electrode. Thus we have 

(2.3) p^ds = Ii, Z = 1,2,...,L, r 

where Ii is the current sent to the ith electrode and ei denotes the part of d that 
corresponds to the Zth electrode. This is combined with 

du 

(2.4) — =0 in the gaps between electrodes. 

The conventional way to model the very high conductivity of the electrodes is 
to impose the constraint that u is constant on each one. These constants, which we 
denote by Vi, are the voltages we measure. We write these constraints as 

(2.5) u="Vl"onei, ' / = 1,2,...,Z, 

where the quotes are used to remind us that the VJ are not speci ed in advance but 
are part of the solution of the forward problem. This model we call the shunt model. 

The shunt model, unfortunately, doesn't give results that agree with experimental 
data either [CING]. It fails to account for an electrochemical e ect that takes place 
at the contact between the electrode and the body. This e ect is the formation of a 
thin, highly resistive layer between the electrode and the body. The impedance of this 
layer is characterized by a number zi, which we call the e ective contaict impedance 
or surface impedance. We therefore replace constraint (2.5) by 

(2.6) u^zi ^ = wr on ei, / = 1, 2, 3, . . . , L. 

The resulting model we call the complete model. 

The complete model consists of (2.1), (2.2), (2.3), (2.4), and (2.6), together with 
the conditions 

L 

(2.7) ^ Ii = 0 (conservation of charge) 
i=i 

and 

L 

(2.8) ^ Vi =0 (choice of a ground). 
i=i 

This model has been shown to have a unique solution [SCI]. It is able to predict the 
experimental measurements to better than .1%. 
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3. Some Issues of System Design. At present, there is only one commercial 
EIT system. This is the applied potential tomography (APT) system invented by 
David Barber and Brian Brown [BB]. This system uses a single current source and 16 
electrodes to make low-resolution images of conductivity changes inside the body. 

One way to improve this system would be to improve its resolution. How can this 
be done? Certainly better resolution must involve using more electrodes. To study 
this idea, we introduce a quantitative measure of the ability of a current density 
j to distinguish two di erent admittivities and . Distinguishing between them 
means that the voltage di erence \\R{ )j R{ is greater than our measurement 

precision. (Here || ||oo denotes the supremum.) However, this voltage di erence can be 
made arbitrarily large by applying a su ciently high current. Since it is not practical 
or safe to apply all the power in the universe, we need some constraint. 

What constraint is appropriate? Researchers have variously constrained the max- 
imum magnitude of the currents [CI], the sum of the currents [EP, KE, KVKK], and 
the power [CI]. The present safety regulations [AAMI] constrain the sum of the cur- 
rents. However, if a current of a given magnitude is applied to a su ciently small 
electrode, the current density can become so high as to cause pain. These regulations 
are thus not appropriate, and are therefore being reexamined. 

For medical applications, there are two physiological e ects to consider. At low 
frequencies, nerves and muscles can respond to electric currents. At high frequencies, 
organs are subject to damage by heating. We believe that both safety concerns can 
be addressed by constraining the applied current density, which bounds the applied 
power as well. 

Accordingly, we de ne [I] the "distinguishability'' of from by a current 

density j to be 

(3.1) 0)- |,^.„ 

where R{ )j denotes the electric potential or voltage on d resulting from the appli- 
cation of the current density j to a body containing the axlmittivity distribution . 
Here, for simplicity, we use the L^{d ) norm. A discussion of other norms can be 
found in [CI]. 

It was shown in [GIN2) that as the area of d on which current is applied shrinks to 

zero, the distinguishability also goes to zero. In particular, if one increases the number 
of electrodes while applying current only between a pair of them, the distinguishability 
decreases. 

This suggests that in order to improve resolution by increasing the number of 
electrodes, one should apply current to all the electrodes. It is for this reason that 
the Rensselaer ACT system has as many current generators as electrodes [GINlj. 

A similar result [GIN2] also implies that EIT systems should use large electrodes 
that 11 as much of d as possible. 

For a many-electrode, multiple-current-generator system, the question arises of 
which current density patterns should be used in order to best distinguish between 
two admittivity distributions and 

We say that the current density i is a "best" pattern for distinguishing from if j 
maximizes the distinguishability . Many readers will recognize (3.1) as the Rayleigh- 
Ritz quotient; it achieves its ma^mum when j is the eigenfunction of |ii( ) ^( )| 
corresponding to the largest eigenvalue. 

As a simple example [GINl], consider the problem of distinguishing a homoge- 
neous annulus from a homogeneous disk with the same conductivity. In this case, 
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the eigenfunctions of ) R{ )\ are trigonometric functions and the best current 
density is therefore simply j{ ) = cos( o) for any o- This example can be used 
to determine the size of the smallest object detectable in the center of an otherwise 
homogeneous medium by measurements of a given precision. We study the center 
because objects placed there have the least e ect on boundary measurements. 

In particular, [GINl] calculated the size of the smallest centered cylindrical insula- 
tor (or conductor) whose presence can be detected in a homogeneous cylindrical tank 
with roughly the same diameter and conductivity as a human chest, when probed by 
a 32-electrode system with .1% accuracy. The results of these calculations agree with 
the experiments [CING, GINl]. In particular, they show that a single current applied 
between an adjacent pair of electrodes cannot distinguish objects smaller than st size 
in the center of a chest-sized tank. A single current applied between diametrically 
opposed electrodes cannot detect objects smaller than two ngers in diameter, while 
the "best" cosine pattern can detect an object the diameter of a single nger in t^e 
center of a chest-sized tank. 

If the conductivity is not rotationally invariant, then the cosine is not necessarily 
the best current density to distinguish from a homogeneous conductivity . In 
general, since the best current densities depend on the unknown conductivity inside 
the body, they cannot be known in advance. They can be determined, however, by 
the adaptive process described in Appendix 2. 

The alert reader may be wondering at this point why we pay so much attention to 
best current densities. After all, the mathematical model is linear, so on a system with' 
a limited number of electrodes, we should be able to apply any linearly independent 
set of current patterns: and from the corresponding measured voltages, synthesize the 
result of applying any other current pattern. 

The problem is that the measurement process introduces nonhnearities. For exam- 
ple, a nonzero voltage that is smaller than the measurement precision of the voltmeter 
will register as zero. In other words, the measurement process causes information to 
be lost. An adaptive measurement scheme for obtaining the most possible information 
is discussed in Appendix 2. 

How many current patterns should we apply? Although in theory, we would need 
to apply all possible current patterns in order to obtain all possible information, in 
practice this would take too much time. Because the problem is very close to being 
linear, we apply only a linearly independent set of patterns that are chosen to have 
maximal information content (in the sense of being the best for distinguishing the 
medium from a homogeneous guess). 

4. Reconstruction Algorithms. The reconstruction problem is to obtain an ap- 
proximation to in the interior from the boundary measurements. This problem 
is challenging because it is not only nonlinear, but also ill posed, which means that 
large changes in the interior can correspond to very small changes in the measured 
data. 

From a theoretical point of view, all possible boundary measurements do uniquely 
determine the admittivity in the interior [KV, SU, N]. However, in practice we are 
limited to a nite number of electrodes and a nite number of current patterns. 

Many reconstruction algorithms have been proposed. We outline several of the 
di erent approaches and then describe in more detail the methods we used to make 
the images that accompany this paper. 

The di erent approaches fall into several categories. The rst are based on linear 
approximations. These are noniterative methods based on the assumption that the 
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conductivity does not di er very much from a constant. Examples of linear methods 
are the Barber-Brown backprojection method [BB] and related methods [SV, BT], 
Calderon's approach [C, IC2, 11, CII, IC], moment methods [BAG, CW, AS], and 
one-step Newton methods [CINGS, B, ESIC, FCIGN, G, S]. The images in this paper 
were made with a one-step Newton method that is close to being an approximate 
linearization; this we discuss in more detail below. 

Another class of methods is iterative methods. These include, typically, output 
least squares for various functionals. Examples include [ESIC, KM, WFN, YWT, 
JIEN, H, WWT, BP, J, D, K, SV]. The papers [D] and [K], in particular, contain 
rigorous proofs of convergence of Newton-type methods under various conditions. 

A related class includes the adaptive methods, in which the applied patterns of 
current are adjusted to get the best reconstruction [GINl, GIN2, NGI, S, BP, ICp]. 

One promising method is a layer-stripping algorithm [Sy, SOU]. It is based on the 
idea of rst nding the electrical parameters on the boundary of the body, then math- 
ematically stripping away this outermost known layer. Then the process is Tepeated,: 
and the medium is stripped away, layer- by- layer, with the electrical parameters being 
found in the process. This method is appealing because it is fast, addresses the full 
nonlinear problem, and works well on continuum-model synthetic data. However, no 
available layer-stripping algorithm works well on complete-model data. 

There are also more theoretical papers that present formulas or suggestions for 
reconstructions, implementations of which have not been published. Such papers are 



5- The IMoser Algorithm and Its Implementation. For an L-electrode system, 
we apply a linearly independent set of current patterns. Because of the constraint 
(2,7) , a full set contains L 1 current patterns, which we denote by , . . , , ^ , 
where = (/f,/2, • j^i)- The corresponding voltage patterns we denote by V^, 
y^, . . . , ^. One might think that from voltages measured on L electrodes, for L 1 
current patterns, we would have L{L 1) independent degrees of freedom. However, 
by (2.8), the voltages are constrained by the choice of a reference voltage, so in fact 
we have only L 1 measurements for each current pattern. Moreover, some of this 
data is redundant, because the current- voltage map is symmetric. Thus the number 
of independent degrees of freedom is the number of degrees of freedom of a symmetric 
{L 1) (L 1) matrix, namely, L{L l)/2. For our ACTS system, which has L = 32 
electrodes, the number of degrees of freedom is 496. 

For simplicity, we consider only the problem of reconstructing the conductivity , 
or, alternatively, the resistivity =1/ . We can only hope to recover a limited number 
of degrees of freedom of ; we denote these degrees of freedom by i, 2> • • • » Ni where 
TV L(L l)/2. For example, could be represented by specifying its averages over 
small mesh elements; the average over the nth mesh element would be denoted 
Henceforth we will assume that is completely determined by the „ . We would like 
to nd a conductivity whose voltage patterns ),U'^{ ), - ., ^( ) are equal 
to the measured ones. To do this, we attempt to minimize the functional 



[SU, N, R, R2, C, 102] . 



L 1 



(5.1) 





To minimize this functional, we di erentiate with respect to each degree of freedom 
n, and set the derivatives to zero. This gives us the following set of N nonlinear 
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equations to solve: 

A standard method for solving such a system of nonlinear equations is Newton's 
method; the NOSER code takes one step of a (regularized) Newton's method, from an 
initial guess of a uniform resistivity, and displays the result. The advantage of taking 
only one step is that the Jacobian matrix for the uniform case can be calculated ahead 
of time and stored. Thus the reconstruction is done via the formula (from [E]) 

(1) = °c+2x:e''m5«='^ 

where C and 5*'-^ are precomputed vectors of length 496 that depend only on known 
factors such as the geometry and regularization, and where ° and Vkj are scalars 
computed from inner products of the data with known vectors. 

This algorithm was implemented on a 386 personal computer with an Alacron 
AL860 AT printed circuit board, which relies on the Intel i860 microprocessor. This 
board has a data retrieval rate of 160 Megabytes per second and is capable of 80 
Mega ops maximum throughput. Our implementation of the above algorithm on 
this board achieves a rate of 60 reconstructions. per second. More details about the 
implementation are given in [E]. 

6. Some Experimenta! Tests. Our present third generation system, ACTS, in- 
corporates the fast algorithm of section 5 together with electronics that are designed 
according to the principles of section 2. This system is capable o f ma king 20 images 
per second with data accurate to one part in 2^^ , using a single array ot '62 electrodes. 

We have used this system to conduct many experiments, some of which are illus- 
trated in the gures. We have done experiments in a test tank (Figures 1 and 2), on 
a human volunteer (Figures 3, 4, and 5), and on dogs (Figures 5 and 6). 

The dog experiments are intended to test the feasibility of using electrical imped- 
ance tomography to detect pulmonary emboli. The dog was prepared so that the left 
and right lungs could be ventilated separately. Figure 6 shows the result of closing 
o the ventilation of one lung at a time. Then a balloon catheter was used to occlude 
a major branch of the pulmonary artery, to simulate a pulmonary embolus. Figure 7 
shows the result. 

7. Future Challenges. Many challenges still remain to be overcome before EIT 
will provide a clinically useful device. The challenges can be classi ed into the areas 
of electronics, algorithms, and clinical applications. 

Because the problem is intrinsically so ill posed, the measurements must be made 
very accurately. For clinical applications, they must also be made fast. Finally, a 
great many electrodes should be used to obtain good images. In particular, because 
current does not con ne itself to a plane, measurements should be made on as large 
a surface as possible, which means that many electrodes should be used. 

The ill-posedness and nonlinearity of the reconstruction problem also make algo- 
rithm development di cult. The NOSER and the backprojection algorithms have the 
advantage of being fast, but they are not very accurate. Algorithms need to be devel- 
oped that are fast and accurate and that apply to a wide range of surface geometries. 
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Fig. 2 Images of the resistivity of two di event test tanks like the one shown in Figure 1 The two 
dz event tanks had hearts of di erent sizes, meant to simulate di erent times during the 
heart s cycle. 



Finally, the ill-posedness of the problem also makes it unlikely that EIT images, 
even with many electrodes, will have resolution comparable to that of CT or MRI 
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Fig. 3 The ACT 3 system with. 32 electrodes encircling the chest of a subject. Here the goal is to 
make images that show the changes in volumes of air and blood that occur with breathing 
and the pulsatile circulation of the blood. Such images are referred to as ventilation and 
perfusion images, respectively. This is one possible positioning of electrodes for this purpose. 

images. However, EIT is low cost, noninvasive, and provides information about the 
electrical parameters of the body, which is information that cannot be obtained by 
these other methods. Speci c clinical applications still need to be explored. 

Appendix 1 : Derivation of (2.1) from Maxwell's Equations. The xed-frequency 
version of Maxwell's equations is 



(A.1) 
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Fig. 4 This is a 'perfusion" image of a human subject. I t'waa' produced- vtith the electrodes con gured 
as in Figure 3. This image a di evence image, forrned'by.subtract^g images taken at two 
di erent times. The rst image wtis made when the heart hs ventxicles have contracted, and 
the second was made when the heart's ventricles are lied with blood. In this di erence image, 
the lung region shows an increase in the magnitude of the admittivity, while the ventricles' 
region shows a decrease in admittivity. Thijs is' because blood (which has high admittivity) 
has traveled from the heart to the lungs. 



(A.2) VAH= E-^iivE, 

where E denotes the electric eld, H denotes the magnetic eld, and VA denotes the 
curl operator. In order to determine whether our parameter ranges are such that we 
can nd a simplifying approximation to these equations, we rst write the equations 
in nondimensional form. To this end, we write E == [E]E, H = ^ = [a:]x, and 

VA = [x] ^VA, where the quantities in brackets are scaiars carrying the units, and 
the quantities with tildes are nondimensional vectors. 
With this notation, we can write (A.l) and (A.2) as 

(A.3) VA^= iu 

(A.4) VA^= imE + i^l^E. 

If we now choose units for [E] and [H] so that [-fi^JM/fif] = 1, then we can write 
(A.3) and (A.4) as 

(A.5) VaE= iuj 

(A.6) V A ^ = E -h i—E. 
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Fig. 5 These graphs illustrate the periodic lling and emptying of the heart ( top) and lungs (bottom). 

The top curve is an average of the magnitude of the admittivity over one of the ventricular 
regions shown in Figure 4^ plotted against time. The bottom curve is an average of admittivity 
of the lung region shown in Figure 4, plotted against time. Note that the admittivity of the 
ventricular region decreases rapidly when the ventricles contract and pump highly conductive 
blood out to the body, while simultaneously the lungs^ admittivity increases rapidly as they 
II with blood. 

Electrical impedance tomography systems operate in the range where uj [x]^ is neg- 
ligible. The ACTS system, for example, operates at 28.8 kHz, and is applied to bodies 
smaller than 1 meter in which the conductivity is generally less than 1 (Ohm-meter) ^ . 
The system of [ESIC] , although used to examine metals of high conductivity, operated 
at a very low frequency. 

If the right side of (A.5) is negligible, then we can neglect the right side of (A.l) 
to conclude that E is the gradient of a potential. In particular, we write 

(A.7) E = Vu, 

where u is the electric potential. Using (A.7) in the equation obtained by taking the 
divergence of (A. 2), we obtain (2.1). 
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Fig. 6 Mean admittivity magnitude in the two regions of interest shown in the diagram at the top 
is shown versus time in the three panels at the bottom. In the rst panel, a breath with a 
tidal volume o/ 300 ml/side was iidministered to BOTH lungs: In the middle panel, a 300 ml 
breath wets applied to only the LEFT lung; in the third panel, the breath was applied to only 
the RIGHT lung. The curves in each panel show the admittivity magnitude in the region of 
the left chest at the top of each panel, and in the right chest at the bottom of each panel. 
The area of the left and right regions of interest was 7.9% and 12.1% of the total image 
area, respectively. The boundaries of these r^egions are iso- admittivity contours selected in 
an image of bilateral ventilation. 



To obtain (2.2), we must consider also a current applied to the surface of the 
body. Thus (A. 2) is modi ed as 

(A.8) VAii' = J^PP'+ E. 

Taking the divergence of (A.8) and using (A. 7) gives us 

(A. 9) V Vw = V J'^PP^. 

We integrate both sides of (A.9) over a pillbox of thickness enclosing part of the 
boundary, and use the divergence theorem. In the limit as the thickness goes to 
zero, we nd 

(A. 10) outdu i^du= Jfr', 

where denotes the outer unit normal vector. We assume that outside the body is 
negligible, so the rst term on the left side of (A. 10) vanishes. We also assume that 
the current J^^w' is applied only on the outer surface of the body, so that J^^^ is 



r 
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Fig. 7 Change in mean admittivity magnitude in two regions of interest with de ation of a balloon 
located in a branch of the pulmonary artery. The regions of interest are shown at top, and 
are contours of iso-admittivity. The reference admittivity of each curve is di erent, and has 
been subtracted to allow both curves to be displayed on the same axes. 



zero. We de ne j = jappi^ ^YiqX j denotes the applied current entering the 
body. This gives us (2.2). 

Careful computations [Doer] show that for the operating parameters of the Rens- 
seleier ACTS system, modeling error is about fteen hundredths of a percent. The 
electronics of ACTS are accurate to about three hundredths of a percent. This is 
accurate enough to detect the modeling error, which implies that the images might 
be improved by using the full Maxwell's equations. 

Appendix 2: Adaptive Process for Finding Best Patterns. In order to nd 
the best single-current density for distinguishing from , we can use the following 
adaptive process [Ip, IC]. 

1) Guess any current density jo for which jo = 0 and ||joll = 1- Set k = 0, 

2) Measure the voltage on 



= R( )3k- 
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3) Compute the voltage on d : 

= )3k. 

4) Compute the new estimate jk-\-i to the best current density by 

jfc+i = ~ — — 



5) If the change in j is less than the measurement precision , i.e., 

Ib'fe+i 3k\\ <,..... 

stop; otherwise increment k and repeat, starting with step 2. 

This algorithm is essentially the power method for nding the largest eigenvalue 
and corresponding eigenvector of a matrix [IK]. 

Numerical and experimental tests of this adaptive process by 32-electrbde ACT. 
systems [NGI, GINl, ESICJ have shown as miich as a 30- fold improvement in distin- 
guishability within 5 iterations. 

Now, suppose we wish not merely to determine whether is distinguishable from 
, but to obtain all possible information that will enable us to distinguish between 
the two, perhaps for the purpose of forming an image of the di erence . We 

could do this by using the above algorithm to nd the best pattern, then searching 
the orthogonal complement for the next best pattern, etc. Alternatively, we can nd - 
all the best patterns at once by the following procedure. 

We denote by §R the boundary map for the complete model, i.e., 9^/ ~ V.^ where / 
and V are now L-dimensional vectors of currents and voltages. We denote by the 
di erence map ?R( ) ^{ )- We want to determine the eigenvectors of without . 
knowing in advance. 

Suppose we rst apply any orthonormal set (T^, Z = 1,2,...,L l}of current 
patterns, and denote the corresponding voltage patterns by {( 3fi)T', i = 1, 2, . . . , Zr 
!}• 

We wish to nd vectors / such that 
(A2.1) ( = 

where the eigenvalue is a scalar. We assume that we can write / as a linear combi- 
nation of the T's: 

L 1 



(A2.2) I=Y,qiT'; 



we need only determine the g's. On both sides of (A2.1) we use the expression (A2.2) 
for /, thus obtaining 

(A2.3) Y.qi{ 5ft)T'= xig.T'. 

/=1 Z=l 

We now take the inner product of (A2.3) with the vector T^, using the orthonormality 
of the T's. This results in 

L 1 

(A2.4) 3i)T')= gfe, 

1=1 
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where we have used ( , ) to denote the inner product. Equation (A2.4) tells us that 
the desired q^s are the eigenvectors of the matrix whose elements are (T^, ( 3f^)T'); 
this matrix, moreover, is one that we can compute from experimental measurements! 
The subtle point is that this process must be iterated. This is because the operator 
is not known exactly; our knowledge of it depends on the current patterns we apply. 
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